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A RELATION OF THE EUROPEAN DATUM TO A 
GEOCENTRIC REFERENCE SYSTEM 

J. G. Marsh 
B. C. Douglas 
S. M. Klosko 


ABSTRACT 


Over 31,000 precision reduced optical observations of GEOS-I 
and II in 70 two-day orbital arcs have been used at Goddard Space 
Flight Center (GSFC) in a dynamical solution to determine center- 
of-mass coordinates for 15 tracking stations on the European Datum. 
Comparisons with the results obtained at Centre National d'Etudes 
Spatiales (CNES) give agreement of about 1.5 ppm for chord lengths. 
After considering a scale correction to the European Datum (ED) of 
1950 to account for the absence of geoid heights at the time of its 
reduction, agreement to a few ppm between the CNES/GSFC and the 
ED chords is obtained. However, a small systematic difference 
between survey and satellite results remains for stations in south- 
eastern France and Switzerland. 
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THE RELATION OF THE EUROPEAN DATUM TO A 
GEOCENTRIC REFERENCE SYSTEM 


I 

INTRODUCTION 

Prior to the work presented here, two investigations of geodetic parameters 
on the European Datum using optical flash data had been published. Gaposchkin 
and Lambeck (1970) at the Smithsonian Astrophysical Observatory (SAO) per- 
formed what was essentially a geometric solution for optical stations and two 
French lasers in Europe having first dynamically estimated certain Baker-Nunn 
camera positions in the area. Lambeck (1971) later used these results to deter- 
mine the orientation and scale of the European Datum with respect to the global 
center-of-mass system. At Centre National d’Etudes Spatiales (CNES) in France, 
Cazenave and her associates (1971) used purely geometric techniques with opti- 
cal and laser data to recover the chord distances between San Fernando, Spain 
and eight other sites in Europe. By comparing their chord lengths with those of 
the surveys, CNES determined a scale parameter for the European Datum. The 
scale recovered by CNES differs from that of Lambeck (1971) by 20 x 10" 6 . The 
results described here agree with those of CNES to better than 1. 5 x 10' 6 . 

In our work, the center-of-mass coordinates of the European optical sites 
have been estimated dynamically. In contrast, the French sought relative sta- 
tion positions, while Lambeck’s use of center-of-mass Baker-Nunn positions 
held constrained in his geometric recovery permitted the entire European Datum 
to be estimated in the geocentric system. The present solution used a total of 
70 two-day arcs of GEOS-I and II optical data in a simultaneous dynamical solu- 
tion for the station coordinates . 

All of the European stations except Malvern and Winkfield, England were 
allowed to adjust independently in our dynamical solution. Due to their proxim- 
ity and the small amount of Malvern data available , Malvern and Winkfield were 
constrained to adjust in parallel. The arc length of two days gives enough data 
(an average of 480 observations/arc) for the dynamical determination of the 
satellite orbit without gravity model errors becoming excessive. The arcs were 
selected to optimize the tracking geometry, with those stations tracking both 
GEOS-I and II having data on all sides of the station and in opposing directions. 
This- gave a beneficial cancelling of model error effects. For the stations which 
tracked GEOS-II alone, this requirement could not be met because of satellite 
viewing conditions. Three conditions helped reduce the effect of model errors 
for these GEOS-II stations. First, the strong presence of the well-determined 
GEOS-I and GEOS-II stations at San Fernando, Haute Provence, Winkfield, 
Malvern, Naini Tal, and Addis Ababa in the solutions kept the orbital error 
small over Europe. Second, the arcs were selected so that the well- determined 
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stations in the Republic of South Africa (UOBUR and lOLFAN) and Madagascar 
(1TANAN) were tracking on the same satellite revolution as that of the European 
optical stations. This selection of arcs containing South African optical data 
helped to reduce the satellite position error over Europe. As important as 
either of these steps is that the GEOS-II stations had large amounts of data si- 
multaneous with other stations for which good coverage existed. 

Table 1 presents the number of observations in our final solution for the 
European stations which are shown in Figure 1. In addition, about 25, 000 ob- 
servations from a world-wide network of stations (held fixed) (Marsh, Douglas 
and Klosko, 1971) were used in the two day arcs. The values obtained in this 
solution are presented in Table 2. Their corresponding survey positions on the 
European Datum are presented in Table 3. 


EVALUATION OF THE RESULTS 

Table 4 presents the translation, scale, orientation parameters and correla- 
tion coefficients relating our dynamical positions to those on the uncorrected* 
European Datum. Table 5 presents a list of the stations used with associated 
residuals for this solution. As readily seen in Table 4, all of the rotation pa- 
rameters are very highly correlated with the translation parameters. The 5.2 
x 10 -6 value for the scale for the uncorrected datum obtained by GSFC agrees 
very well with the CNES value, the difference of approximately 2. 5 x 10~ 6 being 
primarily due to the Riga and Greece chords (Table 6). These results are in 
disagreement with the scale value of -12.4 x 10~ 6 obtained by Lambeck from the 
SAO solution. Comparison by Vincent, Strange, and Marsh (1971) of SAO sta- 
tion heights with the SAO (1969) Standard Earth geoid indicated a systematic 
height difference of about 20 meters for the European stations. This difference 
and the scale result may be related. 

Table 6 compares the chord lengths obtained from the GSFC dynamical solu- 
tion with the geometric solutions of SAO and CNES. The agreement between 
GSFC and CNES is very good, 6 of 8 chords agreeing to better than 3. 5m. This 
result is especially significant because of the very different estimation techniques 
employed by CNES and ourselves. The disagreement of the chords to Riga, 
Latvia, and Dionysos, Greece deserves mention. In the case of Riga, Latvia, 
the same local survey was used by GSFC and CNES, but the accuracy of the sur- 
vey was not available. Therefore, the survey value cannot be used to resolve 
the GSFC-CNES disagreement. The disagreement at Dionysos, Greece is prob- 
ably caused by Greece being on the periphery of the geometric net and therefore 


*Uncorrected for scale introduced by absence of geoid height information (Bomford 1971). 
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being constrained in only limited directions in the CNES solution; i.e. , a ma- 
jority of the usable data is on one side of the station. However, the mean agree- 
ment between the CNES and GSFC results , including these two questionable 
sites, is still only five meters . 

Due to scarcity of data, the GSFC solutions for Oslo (Norway), Naini Tal 
(India) and Meudon (France) are weaker than the others. The GSFC solution in 
particular, differs by about 40 m in the Z component from the Naini Tal solution 
given by SAO. Further investigation is required to resolve this discrepancy. 
Delft (Netherlands) also had limited data in our solutions, but the geometry is 
excellent and the data is simultaneous with other European stations. The uncer- 
tainties in the ties to the European Datum for Riga, Latvia and Addis Ababa, 
Ethiopia, are reflected in the results. These sites are believed to have strong 
dynamical solutions. 

According to Bomford (1971), the European Datum contains a systematic 
scale error due to the unavailability of the geoid heights throughout this system 
at the time of its reduction in 1950. When the baseline distances were reduced 
to the International Ellipsoid, the geoidal height variation in Europe from 
Potsdam was of necessity omitted. With more recent work on the European 
continent and the availability of geoid heights, distances in Europe should be 
corrected by +lppm for every +6. 4m by which the geoid is actually above the 
spheroid. When the GSFC and CNES chords from San Fernando are compared 
with survey chords corrected in this manner, the large discrepancy between 
satellite and surveyed distances is reduced. These corrected values are shown 
in Table 6. 

As seen in Table 6, the chord lengths to San Fernando obtained from the 
SAO Standard Earth (Gaposchkin and Lambeck, 1970) are in disagreement with 
those obtained by CNES and GSFC. The SAO chords are smaller than the sur- 
veys, the satellite-determined chords are larger for the GSFC and CNES solu- 
tions . While the correction indicated by Bomford removes much of the disagree- 
ment between CNES and GSFC chord lengths compared to those of the surveys, 
the SAO scale factor becomes further in disagreement. Therefore, no further 
comparisons are made in this text between SAO chords and those of the corrected 
surveys. 

Table 7 presents a comparison of the differences between the chord lengths 
of GSFC with corrected and uncorrected ground surveys throughout the European 
Datum. Those chords which are clearly doubtful are shaded;. DEZEIT (Addis 
Ababa, Ethiopia) and RIGALA (Riga, Latvia) are tied to the European Datum 
with uncertain accuracy. As shown in Table 7, surveyed chords corrected for 
scale as recommended by Bomford are in better agreement with the GSFC satel- 
lite solution. Table 8 presents a similar comparison for SAO considering only 
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the uncorrected ground survey. Figure 2 presents a series of histograms of 
GSFC and SAO's chords compared with those of the European Datum, also show- 
ing GSFC compared with the corrected surveyed chords . The GSFC chords 
agree with the corrected ED chord lengths to better than 10 meters in 44 of 66 
cases. This is consistent with our accuracy estimates for our satellite solutions 
which indicated that our recovered station locations are normally accurate from 
between 2 to 7 meters in each coordinate. 

Concerning Shiraz, Iran, and Naini Tal, India, according to Bomford (1971) 
these sites were tied to the European Datum in the middle 1960’s with geoid 
height information available. Therefore, chords to these two sites require no 
additional scale correction related to geoid height. 

Figures 3, 4 and 5 present the shifts in the rectangular coordinates for the 
dynamically recovered stations from the uncorrected European Datum to the 
geocentric system. Note that the Baker-Nunn camera position at San Fernando, 
Spain, is systematically different in AX from those sites nearest to it. We 
further note that the station positions located in southeastern France and Switzer- 
land seem somewhat inconsistent in AX, with systematic differences of about 
10m apparent. This systematic variation in AX in central Western Europe is 
not explained by the geoid height offset in the area. Due to a scarcity of station 
positions, this result is not conclusively demonstrated throughout Europe along 
this parallel. 


CONCLUSIONS 

This paper demonstrates that satellite geodesy can approach the level of 
accuracy long associated with classical surveying techniques. Other investiga- 
tors in the past have used satellite data to connect isolated tracking stations 
with major geodetic systems , but this is believed to be the first time that satel- 
lite solutions have successfully been used to detect systematic errors within a 
major geodetic datum. Dynamic satellite techniques provide the advantage that 
large surface areas can be adjusted simultaneously with an accuracy almost in- 
dependent of distances between stations . Previous solutions by CNES and SAO 
did not consider a scale correction to the survey chords due to the absence of 
geoid height information at the time of the reduction of the European Datum in 
1950. Agreement of the GSFC and CNES solutions with the survey improved 
significantly when this correction was added. A systematic difference remains 
for stations in southeastern France and Switzerland. Comparisons of GSFC's 
results with those of CNES give general agreement of 1. 5 ppm or better for chord 
lengths and suggests that the previously published SAO solution, which differs 
by 20 ppm in scale from the French results, contains systematic errors. When 
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surveyed lengths are corrected as suggested by Bomford, the GSFC solution 
agrees with surveyed distances to a few ppm. 
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Table 1 


Number of Optical Observations per Station 
Used in Dynamical Solution 


LOCATION 

CODE 

NAME 

STATION 

NUMBER 

OBSERVATIONS 

Winkfield, England 

1WNKFL 

1035 

611 

Delft, Netherlands 

DELFTH 

8009 (9065) 

144 

Zimmerwald, Switzerland 

ZIMWLD 

8010 (9066) 

481 

Malvern, England 

MALVRN 

8011 (9080) 

87* 

Haute Provence, France 

HAUTE P 

8015 

779 

Nice, France 

NICEFR 

8019 

999 

Meudon, France 

MUDONI 

8030 

203 

San Fernando, Spain 

1SPAIN 

9004 

1750 

Naini Tal, India 

1NATOL 

9006 

161 

Shiraz, Iran 

1SHRAZ 

9008 

41* * 

Addis Ababa, Ethiopia 


9028 

337 

Dionysos, Greece 

GREECE 

9091 

1027 

Oslo, Norway 

OSLONR 

9426 

28** 

Riga, Latvia 

RIGALA 

9431 (9074) 

453 

Uzhgorod, U.S.S.R. 

UZHGOR 

9432 (9077) 

395 


*1MALVRN was held constrained to 1WNKFL 

**Only one right ascension and declination observation was precisely reduced for each pass of data from 
these stations in this time period. 
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Table 2-a 


Estimated Station Coordinates (X, Y, Z) 


Station 

X 

Y 

Z 

Name 

Number 

(M) 

(M) 

(M) 

1WNKFL 

1035 

3983102 

-48512 

4964720 

DELFTH 

8009 

3923391 

299885 

5002982 

ZIMWLD 

8010 

4331307 

567522 

4633122 

MALVRN 

8011 

3920151 

-134739 

5012737 

HAUTE P 

8015 

4578335 

457982 

4403200 

NICEFR * 

8019 

4579471 

586614 

4386422 

MUDONI 

8030 

4205620 

163727 

4776555 

1SPAIN 

9004 

5105586 

-555238 

3769681 

1NATOL 

9006 

1018208 

5471117 

3109585 

1SHRAZ 

9008 

3376880 

4403985 

3136261 

DEZEIT 

9028 

4903769 

3965210 

963853 

GREECE 

9091 

4595174 

2039458 

3912663 

OSLONR 

9426 

3121268 

592634 

5512724 

RIGA LA 

9431 

3183873 

1421477 

5322789 

UZHGOR 

9432 

3907419 

1602436 

4763906 















































































Table 2-b 


Estimated Station Coordinates (0, X, h)* 


Station 

Geodetic Latitude 

East Longitude 

Ellipsoid 

Height 

(Meters) 

Name 

Number 

Deg 

Mn 

Second 

Deg 

Mn 

Second 

1WNKFL 

1035 

51 


46.40 

mu 

18 

7.93 

90 

DELFTH 

8009 

52 

H 

6.76 

E 

22 

15.29 

46 

ZIMWLD 

8010 

46 

52 

37.18 

B 

27 

| 

933 

MALVRN 

8011 

52 

8 

36.42 

358 


na 

137 

HAUTEP 


43 

55 

57.55 

5 

42 

44.74 

694 

NICEFR 


43 

43 


B 

17 

58.5*8 

405 

MUDONI 

8030 

48 

48 


2 

13 

45.94 

190 

ISP AIN 

9004 

36 

27 

46.99 

353 

47 

36.31 

55 

1NATOL 

9006 

29 

21 

33.31 

79 

27 

27.07 

1856 

1SHRAZ 

9008 

29 

38 

13.80 

52 

31 

11.25 

1564 

DEZEIT 

9028 

8 

44 

50.71 

38 

57 

32.98 

1901 

GREECE 

9091 

38 

4 

44.39 

23 

55 

58.43 

490 

OSLONR 

9426 

60 

12 

39.50 

10 

45 

2.69 

595 

RIGALA 

9431 

56 

56 

55.32 

24 

3 

32.17 

-15 

UZHGOR 

9432 

48 

38 

1.46 

22 

17 

54.88 

205 


*a e = 6378155 meter., f = 1/298.255 
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Table 3 


Station Coordinates on the European Datum 


Station 

X 

Y 

Z 

Name 

Number 

(M) 

(M) 

(M) 

1WNKFL 

1035 

3983202.6 

-48394. 5 

4964835.4 

DELFTH 

8009 

3923486.0 

300006.0 

5003095.8 

ZIMWLD 

8010 

4331390.6 

567637.4 

4633235.9 

MALVRN 

8011 

3920250. C 

-134624. 4 

5012852.2 

HAUTE P 

8015 

4578413. 0 

458091.0 

4403312.0 

NICEFR 

8019 

4579554.2 

586729. 1 

4386535.6 

MUDONI 

8030 

4205717.7 

163840. 9 

4776060.8 

1SPAIN 

9004 

5105680.1 

-555102.9 

3769799.3 

1NATOL 

9006 

1018274.4 

5471244. 5 

3109773.8 

1SHRAZ 

9008 

3376966.7 

4404122.1 

3136407.9 

DEZEIT 

9028 

4903853.4 

3965302.9 

964020.8 

GREECE 

9091 


2039577.4 

3912795.2 

OSLONR 

9426 

3121372.8 

592748. 1 


RIGALA 

9431 

3183998.7 

1421638.2 

5322894.3 

UZHGOR 

9432 

3907494. 2 

1602533.2 

4764034. 8 
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Table 4 


U 

H 



10 


>mega 










































Table 5 


Stations in Transformation Solution 
with Their Associated Residuals 



Residual (meters) 

X 

i 

Y 

Z 

Malvern, England 

3.4 

-2.2 

3.3 

Nice, France 

-3.9 

-5.9 

-4.6 

Zimmerwald, Switzerland 

warn 

-2.4 

mm 

Dionysos, Greece 

1.8 

3.3 

3.4 

Delft, Netherlands 

2.0 

6.7 

-0.6 

San Fernando, Spain 

1.4 

0. 6 

2.3 

RMS of fit 

3.1 

4. 1 

3.3 
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Comparison of Chord Distances from Station 9004 (San Fernando, Spain) 
on the Corrected and Uncorrected European Datum of 1950 
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Table 7 


HAUTEP 

NICEFR 

ZIMWLD 

MUDONI 

DELFTH 

MALVRN 

OSLONR 

GREECE 

UZHGOR 

ISHRAZ 

INATOL 

RIGALA 

DEZEIT 


HAUTEP 

NICEFR 

ZIMWLD 

MUDONI 

DELFTH 

MALVRN 

OSLONR 

GREECE 

UZHGOR 

ISHRAZ 

INATOL 


Differences Between GSFC Satellite and Ground Survey 
Chord Lengths (Meters) 
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RIGALA 



Differences Between SAO Satellite and Ground Survey Chord Lengths (Meters) 
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Stations on the European Datum 
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Figure 2. Histogram of the Difference Between Surveyed and Satellite Derived Chord Lengths 
on the European Datum 
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Figure 4. Adjustments in the Y-Coordinate (Satellite-Survey Solution) nv Meters 



